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ABSTRACT 

The central engine that drives gamma ray burst (GRB) explosions may derive from the 
ability of electrons/positrons and nucleons to tap into the momentum and energy from 
the large neutrino luminosity emitted by an accretion disk surrounding a black hole. 
This transfer of momentum and energy occurs due to neutrino absorption, scattering, 
and annihilation and the non-spherical geometry of the source both increases the 
annihilation efficiency and, close to the black hole, directs the momentum transfer 
towards the disk axis. We present annihilation efficiencies and the momentum/energy 
transfers for a number of accretion disk models and compute the critical densities of 
infalling material below which the transfer of neutrino momentum/energy will lead 
to an explosion. Models in which the neutrinos and antincutrinos become trapped 
within the disk have noticeably different momentum and energy deposition structure 
compared to thin disk models that may lead to significant differences in the explosion 
dynamics. 
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1 INTRODUCTION 

Gamma ray bursts were first observed forty years ago, but only recently has significant progress been made in under- 
standing their origin. Evidence is mounting that at least t he 'long durat i on' bursts are a s sociated with a rare type of 
supernova event such as a 'failed sup ernova' or 'coUaps ar ' JWooslevlHooi |Paczvnskilll99d iMacFadven fc Wooslevlll999l : 
iM^^F ^c^v cn. Wooslev. and Hegeill20 01: Proea et al.' with neutro n-star mergers as a candidate for the shorter bursts 

l)Paczvriski..l991. : .Ruffert fc .TankJl999- Rosswog fc R,amirez-R,uiJ20n2l) . In either case, it is likely that an accretion disk sur- 
rounding a black hole forms that cannot cool efficiently by photon emission because of the high densities and temperatures. 
But with their much smaller cross sections the neutrinos, if any are produced, may escape and cool the disk. 

The electron neutrinos and antineutrinos emitted from the disk are produced from electron capture by protons or positron 
capture by neutrons respectively. Both reactions modify the make-up of the disk with electron capture initially predominating 
as material moves inward toward the black hole. The frequency of these reactions (and their inverse) are strong functions of 
the thermal properties of the disk and grow considerably as material moves towards the black hole where the temperatures 
and densities are highest. In turn, the density and temperature are functions of the mass, Mbh, and spin, a, of the black 
hole, the viscosity, a, of the disk and the accretion rate, rh since these four parameters determine the potential and the rate 
at which material flows radially inward. In particular, as either rh or a increase so does the density of the disk (particularly 
in the hottest region close to the black hole) leading to greater neutrino luminosity. 

The neutrino scattering length within the disk is similarly affected by the thermal properties. If sufficiently high tem- 
peratures are reached the neutrinos and/or antineutrinos become trapped and, consequently, their spectrum is altered. Due 
to the change in nucleon composition as material ffows inward, neutrinos and antineutrinos will have different optical depths 
and will not become trapped at the same radius. For the lower accretion rate disks, m < 0.1 Mq s~^, the density of the disk 
is never sufRcient to trap the neutrinos and so they immediately escape. But for the higher accretion rates the neutrinos do 
become trapped. Once trapped the neutrinos thermalize and the neutrinos emitted from these regions have temperatures of 
a few MeV, the exact value depends on the model, and are functions of the radius. In addition, neutrino trapping also has 
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the effect of restricting neutrino emission to that portion of the disk beyond a 'neutrinosurface' which, given the cylindrical 
symmetry, are roughly toroid in shape. 

Descriptions of disks with low a ccretion rates are given in IPopham. Wooslev fc Frverl il999l) while models which do 
include neutrino trapping a re given inlPiMatteo. Perna fc Naravanr(l200^ Additional calculations of neutrino emission from 
accretion disks are given in lSetiawan et alJ l)2004h . 

After escaping the disk the neutrinos and antineutrinos do not propagate unimpeded but may interact with material 
above the disk transferring energy and momentum. The extent to which their energy and momentum is tapped is the primary 
interest of this paper since the imparted energy may^ power the burst. The transfer of momentum and energy occurs by a 
number of processes t hat we shall discuss in jQl We t ake calculations of the neutrino flux emitted from every point on the disk 
from previous work in lSurman fc McLauehli Our results are presented in JSl^-nd to investigate the effects of each of 

the four parameters needed to describe the disk we compute the energy /momentum transferred in a number of models. The 
parameters for each are listed in Table Q 



2 NEUTRINO INTERACTIONS 

The neutrinos and antineutrinos that emerge from the accretion disk can interact in a number of different ways with the 
material above the disk as they propagate. These interactions include 



Uc + Ve ^ e'^ + e , (1) 

Vc + n +P, (2) 

Ue+P e'^ +n, (3) 

Ve + e ^ i/i + e' i>e + e ^ K + e , (4) 

+ P ^ l^'e + p' i>e + P ^ K + P , (5) 

+ n ^ 1/'^ + n i>e + n ^ u'^ + n . (6) 
We discuss each in turn. 



2.1 Neutrino Annihilation 

Annihilation of neutrino-antineutrino pairs to electron-positrons, reaction Q, will deposit energy at a rate per unit volume, 
dLvv/dV , that is a function of the (cylindrical) coordinates r and 2, symmetry ensuring there is no polar angle dependence. 
The expression for dL^v/dV is 

^ = y Jdn,dn,J J dE^dE, U^'^^E,)Mr„E,) 
{E^ + Ei, 



Ely Ei) 



{Ey Ea \vv -^u\ cTvv} ■ (7) 



The quantities fv^{ri,, E^) and fu^{ri^, Ej^) are the differential neutrino and antineutrino number fluxes p er unit area that 
emerge from radial coordinates r,j, Vi, on the disk. We take these fluxes from lSurman fc McLauehlinl ||2004^ and we refer the 
reader to that paper for details of their calculation. The two radial positions may be expressed in terms of the angles (Oi,, <j)i,) 
and {dv, 4>v) together with r and z, as shown in Fig.0 so that 

r1 = [r cos(f)u + zt&nOv)^ + sm^ (f)u (8) 



= [r cos(f)i} + zt&nOi))^ + sm^ (f);}. (9) 



The term in curly braces in equation Q is a Lorentz invariant quantity so can be evaluated in the center of mass frame and 
expressed in the terms of the invariant Mandelstam variable s. When in the CM frame of the collision s = AEvEi?; for the 
frame of the disk s = 2EuE,-j{l — cos 7) where 7 is the angle between the v and v trajectories, given by 

C0S7 = sin 6v sin 6i?,cos((j)v ~ (j>v) + cos 9v cos On. (10) 
The annihilation cross section is well known fsee lHerrera fc Hacva nl ll989l. for example) and so 



^ Though so me of the required energy m ay come by electromagnetic extraction of energy from the rotating black hole, e.g. through the 
mechanism of lBlanford fc Znaiekl il977^ 
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(11) 



with Cv = 1/2 + 2sin^ ^h^, Ca = 1/2 and 9w being the Weinberg angle, while ao = 4G% ml / {tt h"') . Finally, the two S{6, <j)) 
terms in equation (0 are corrections to the neutrino and antineutrino differential number fluxes per unit area that are required 
in order to take into account the fact that neutrinos and/or antineutrinos within the disk may be trapped. When trapping 
occurs we may divide the disk into optically thick and opticaUy thin zones, the optically thick region being closer to the 
black hole, with radial boundaries of and aj^ for neutrinos and antineutrinos respectively. The values of and aj^, where 
appropriate, are listed in table Q. In doing so we make use of the approximation that the scattering length in the optically 
thick zone is much smaller than the disk scale height so that the transition from one to the other is discontinuous. If the 
solid angle di} along the line of sight, at angles {0,(j}), intersects the optically thin zone then the volume of the disk along 
the line of sight is larger by sec compared to that along a direction paraUel to the disk axis. But if the zone were optically 
thick then only neutrinos from the material within one scattering length of the disk 'surface' can be 'seen', a distance that 
does not vary with the observation angle. Hence the differential number surface fluxes in equation Q, which are defined by 
'collapsing' the disk's volume flux along a direction parallel to the disk axis, need to take into account this effect. We adopt 
Sv{Ov, <j)v) = 0{a.v — r^,) + Qifv — a^) secO^, where Q is the Heaviside step function, and similarly for the antineutrinos. 

The limits on the angle integrals are functions of both the position (r, z) and the inner and outer radii of the disk, Rin 
and Rout- We shall only consider the case for which r < Rout so that the limits on (j) are and 2n. The limits for 6 are then: 



z tan( 



z tan( 



when r < Ri„, 

when Rin <r < Rout, 



\/ Piout ^ '"^ sin^ — r cos ^. 



(12) 
(13) 



2.2 Neutrino Absorption 

Absorption of neutrinos and antineutrinos by the nucleons, reactions Q and IpJ, will transfer both momentum and energy to 
the electrons, positrons and nucleons. If we select a position (r, z) at which to evaluate the transfer from neutrinos emerging 
from (0, Q) then the momentum transfer occurs along the unit vector k'. These vectors and angles are illustrated in figure 
In terms of i,j,k the unit vectors i',j',k' are 

i' = — cos (/)cos 61 i — sin 0COS 61 j — sin k, (14) 
j' = sin i — cos j, (15) 
k' = — cos (jisin ^ i — sin (/)sin ^ j + cos ^ k (16) 

where i is a radial unit vector in the plane of the disk and k the axial unit vector. Due to the cylindrical symmetry the 
integration over can only leave the i and k components of k' as non-zero. Thus we obtain an expressions for the rate of total 
momentum and energy transfer to all particles, f'"'"'^ W^"'^"^ at any given position (r, z) per proton or neutron of the form 



dn, [-sine cos (j!>i + cose k] I dE^ ^^^^^^-^ S^O^, ct>,) E, al''^''\ (17) 

' 47r 



w^^"'^ = dn^ dE.i^^^^i^s4e,,(b.)E.al%'K (is) 



The cross sections are given by 



...M - (i±M^ J, _ „,„ 



(^)^ 



= » (i^) (!^y ,/i - ,20, 

with qa the axial- vector coupling constan t and A is the ne utron-proton mass difference. In these cross sections Wm,W]^ are 
two weak magnetism corrections given bv lHorowit j (|2Qq2) 

Wm = H-1.1— (21) 

mn 

Wm = 1-7.1 — . (22) 

m„ 

The limits on the angular integrals are the same as m v v annihilation. 
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2.3 Neutrino Scattering 

The scattering of neutrinos and antineutrinos by electrons, protons and neutrons, reactions Q through J^, will also transfer 
momentum and energy. The added complication of course is that, in addition to the incoming neutrino energy and angle 
integrals, one must also integrate over the outgoing angles. However due to the cylindrical symmetry any i' and j' terms 
must vanish leaving only the k' component. Again, due to the cylindrical symmetry, there can be no momentum flow in the 
j direction. This leads to momentum and energy transfer expressions due to scattering of 

dn, [-sin6'cos(?ii + cos6'k] /d£;,^^i^%i-^ dQ,',p^,-^, (23) 

/ 47r / dill, 



W = I dn„ I dE„^^^^^^^ f dnl T (24) 



with 

E,.{E,, + M) (l-cos6>:.) 

" M + E,{l-cosei) ' ^ ' 

^ gg(i-cose:.) 

M + £^(l-cos6';.) ^ ' 

where M is the mass of the scattered particle and 8l is the neutrino's scattering angle relative to k'. These expressions show 
that the momentum transfer is always ~ Ei,, whatever the mass of the scattered particle, but that the energy transferred 
becomes very small if M 3> Ei,. 

2.3.1 Neutrino-Electron Scattering 

The cross section for electron-neutrino scattering is again well known f see iBahcalll 1198^ : iHerrera fc Hacva nl ll989l: |Passer3 

I2OO3I for example) and usually expressed in the frame of reference where the electron is initially at rest as 

where Te is the kinetic energy of the recoiling electron. The transformation to da^e/dQ,^, the quantity needed for equations 
1231 1 and 1241 1. is achieved by multiplication of 1271 by {2-k)~^ dTe/d{cos9l) with Te{cos6l) given by equation 1261 . For 
electron-antineutrino scattering Ca — *■ —Ca- We use the approximation that the electrons are initially at rest, however, if 
the temperature of the material becomes of order an MeV, or it acquires large velocities, then corrections due to the finite 
electron momentum would need to be taken into account. 



2.3.2 Neutrino-Nucleon Scattering 

The differential cross sections for these processes are given bv lBurrowj i200ll) 

djir = Te^l^J (^^^j(l+5„cose.), (28) 
dcr,yi, an / E, 



dQ,i, 
where 



5; 



^{—) [{Cv-lf + i9\{CA~lf]{l + 5,cose';) (29) 



l + 3<?i' 

{Cv~lf~g\{CA 



{Cv^iy+igliCA-iy 

The cross sections for antineutrino scattering are identical except for the inconsequential transform of gA ^gA- 



(30) 
(31) 



3 RESULTS 

Table ((2J shows th e neutrino plus antineutrino lumi nosity coming from each of the disk models. Our luminosities are slightly 
smaller than those IPopham. Wooslev fc Frverl l)l999l) for the same models by roughly a factor of 1.6. From examination of the 
table the basic trends with variations of the disk parameters show that the luminosity increases with both the accretion rate 
and spin parameter but that larger black hole masses lead to smaller luminosities. For accretion rates as large as 1 Mq s~^, the 
total luminosity from the accretion disk is as high as lO^'^ergs"^. However, for low accretion rate disks, such as 0.01 Mq s~^. 
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the luminosity in neutrinos plus antineutrinos, L^, coming from the disk can be as low as 10^'erg s~^. The fraction of the rest 
meiss energy radiated from the disk, Lu/rh, also shown in table (|5J, can be as large as ~ 10 per cent when the spin parameter 
and /or accretion rate are large. There is also clearly some interdependance amongst the disk parameters with much larger 
increases seen in Iju of the a = 0.1 thin-disk models when m is increased compared to those models with a = 0.01. Similarly, 
the reduction in luminosity due to increased viscosity is apparently more pronounced when the accretion rate is low compared 
to those cases with moderate accretion rates. The details of the model are also seen to have an effect (models 9 and 10) with 
a doubling of the neutrino luminosity for those that include trapping as compared to the thin disk model. 

Table |21 also shows the total neutrino-antineutrino annihilation luminosity, L^v , integrated over the volume above and 
below the disk. The total power deposited when the accretion rate an d/or the spin parameter is la.rge c an be considerable, 
approaching ~ 5 x 10^^ ergs~^ for m = 10 Mq s~^. In comparison with lPopham. Wooslev fc Frveil il999l) our results for the 
same models are smaller by roughly a factor of 5-6. The difference is mainly due to t he fact that they compute the neutrino- 
antineutrino annihilation rate using Eq. (6-1) in IPopham. Wooslev fc Frveil lll999l) while we use Eq. 0. By comparmg 
the models we notice that the annihilation luminosity varies more rapidly with changes in the disk parameters than 
Increasing the spin parameter from to 0.5 produces a factor of 2 increase in the neutrino luminosity but a factor of ~ 8 
in the annihilation luminosity. Likewise the increase in L^ii in the thin disk models as the accretion rate grows is much more 
pronounced, in some models from factors of ~ 200 for to factors of 10^ for L^i;. The parameter interdependance noted 
above for L^, is again seen in Li^^- Interestingly the inclusion of trapping into the model (models 9 and 10) seems to reduce the 
annihilation luminosity slightly. The efficiency with which the neutrino energy is tapped, L^v/Ljj, is also shown in the table. 
Again this displays a greater variation than does Lv reaching ~ 0.3 per cent for large accretion rates and spin parameters but 
can be as small as ~ 4 x 10~^ if the disk is viscous, the accretion rate is small, or the black hole mass is large. 

Focusing up on model 10. we sho w in figure lEJ the V — V annihilation power as a function of r and z. The result is very 
similar to that o f iRuffert et al.l l)l997l) for neutron star mergers. The bulk of neutrino annihilation occurs in an ellipsoidal region 
extending upwards to z ~ 40 km and outwards to r ~ 90 km peaked in the vicinity of the strongest neutrino-antineutrino 
flux namely the innermost regions of the disk. The oblateness is the product of both the diminishing flux as one moves away 
from the brightest zones on the disk and the cross section's depe ndence upon the ang le between the neutrino trajectories. 
Along the z axis dL^j^/dV varies as dL^j^/dV oc z/{z1 -\~ z^y^'^ fc.f. iRuffert et all (Il997l) with a constant equal to ~ 40 km 
for this particular model^. Along the disk axis the peak power deposition occurs at z*/\/6 ~ 16 km but note that, according 
to the figure, the power deposited at this point is more than an order of magnitude smaller than that immediately above the 
innermost portion of the disk. This fact has also been noted bv lCeorge et alJ (I2OO3I) who used a thin disk model taking into 
account general relativistic effects. The electron positron pairs will thermalize locally so one should expect a large temperature 
and pressure difference between the 'hot spots' immediately around the innermost portion of the disk and the z axis. 

In contrast, the annihilation power along the z axis for the thin disk models, 1-9, varies as dL„ii/dV oc 1/(2* + z^)"^ and 
in model 11 this changes to dLun/dV oc z^ /{z1 + z^^. The difference between models 10 and 11 lie in the contribution of 
the thick and thin portion of the disk to the total antineutrino flux at any given point above the disk: ap = 33.6 km in model 
10 vis-a-vis aj^ = 158.1 km in model 11. At large distances, z 2> z^, all the models seem to follow a similar scaling with z, it is 
the structure close to the black hole that is different. Neutrino opacity not only affects the structure of the disk and the total 
annihilation power but it also manifest itself in where this power is deposited and, consequently, the expected temperature 
and pressure gradients produced upon local electron-positron thermalization will be altered. 

Table|21shows the maximum momentum transfer along the z axis for each of the eleven models due to proton and neutron 
absorption and then neutron, proton and electron scattering. We give the momentum along the z-axis because this is the 
direction in which the jet would emerge but we should note that the momentum transfer at other points above the disk is 
larger than along the axis, see Fig. @. The peak momentum transfer varies as strongly with the disk parameters as does the 
neutrino luminosity and neutrino-antineutrino annihilation power in table |5| A comparison between the entries within each 
model shows that the momentum transferred due to neutrino absorption is larger than that in the scattering reactions but if 
the material along the z axis is primarily leptonic then the scattering of neutrinos by electron/positron will dominate. Across 
different models we see increases in Fk as m grows except for very large m, models 10 and 11, where the trend is reversed 
for the neutrinos; larger momentum transfer with increases in the spin parameter and inverse correlation with the black hole 
mass. The effects of the viscosity are mixed with smaller transfers when the accretion rate is small but larger differences at 
moderate rates. Whether trapping is included also seems to play a role with the momentum transfer along the disk axis being 
larger by a factor of ~ 3 — 4. 

In Fig. (EJ we show the momentum transfer from electron-neutrino scattering as a function of r and z again for our model 
10. At large distances from the black hole the momentum transfer is purely radial but as we approach the black hole the 
non-spherical nature of the source becomes apparent. The largest momentum transfer occurs in the vicinity of the inner region 

A least square fit of our results to the equation dLvv/dV = (Co + Ci z"^)/{z^ + z^)" with m, 71,2;*, Co and Co all adjustable and, 
following the numerical methods used in the computation of equation j2j, assignation of an 'error' proportional to the computed values, 
gives {m, n, z*} = {0.87, 3.77, 42.6 km} with the peak at 14.9 km. 
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of the disk where the neutrino fluxes are largest and there is an oblate, spheroidal region, extending up to z ~ 36 km for this 
model, in which the momentum transfer is directed towards the disk axis. The momentum transferred in the z direction also 
reaches its maximum value (listed in tabled at this point which is, therefore, further out than the peak in dL^j^/dV. 

A rough fit to the momentum transferred due to electron-neutrino scattering along the disk axis^ in model 10 is 
Fk oc /(z't + z^)'^ , whereas antineutrino-electron scattering is closer to^ Fk oc z/{z'^ + z'^)^^^. Again the difference is due to 
the rather small extent of the trapped antineutrino zone. Both neutrino-electron and antineutrino-electron scattering in all 
the thin disk models follow the z/{z'^ + z^)^^^ pattern but in model 11 it is z^ /(z^ + z^)^ that is the better match. Thus 
trapping, and its extent, also appear to alter the details of where, how much, and in which direction electrons gain momentum 
especially in the vicinity of the black hole. 

Finally in Fig. Q we show the momentum transfer from proton-neutrino scattering as a function of r and z again for 
our model 10. The same basic structure is seen and, again, the momentum transferred in the z direction reaches its maximum 
value at z ~ 36 km. Ifowever note that the momentum transferred is somewhat larger than that in e — i/ scattering seen in 
figure The reason lies in the cross sections for each: for neutrino-electron scattering a varies linearly with Eu/m^ but 
in proton-neutrino scattering we have a oc E^/m^. For neutrino energies of, say, ~ 5 MeV the number of proton-neutrino 
scatterings is larger by, roughly, an order of magnitude. As we mentioned above, the momentum transfer per scattering for each 
process is comparable so the net momentum transfer, the scattering rate multiplied by the average momentum transferred, 
is greater for proton-neutrino than for electron-neutrino scattering. This would change if the electrons had kinetic energy or 
temperature of order ~ MeV also. 



3.1 Critical Densities 

The energy deposited into electron/positron pairs a nd the s cattering o f neut rinos will retard the inward motion of material 
towards the black hole. Following the arguments in lFrver"fc Mcszarc3 1I2OO3I) for the case of material falling along the disk 
axis, as a mass element, with density p, falls towards the black hole it experiences an acceleration due to gravity, Fq, and 
a deceleration from both scattering plus absorption, Fv, and from neutrino annihilation, F,jv. If the mass elements freefalls 
from infinity then after it has fallen to the height z the change in kinetic energy is 

^v^{z)= / dz [Fa + F, + F,i,) (32) 

The force of gravity is simply Fq — —Gn Mbh p/z^ , we take the force due to neutrino-antineutrino annihilation as Fjyn = 
1/c dLvn/dV , and the force due to scattering and absorption is 

F^ = n„ (F^"'"^ + F„, + F„^) + Up (F^"*"*' + Fp, + Fp^) -I- (F,, + F,^) (33) 

where rn are the number densities. Using Ue — rip, p — {n„ +np) m,u, where rriu is the atomic mass unit, and introducing the 
neutron fraction Y = nn/ijin +np), then 

F. = fy F„ + (1 - Y) {Fp + F,)l (34) 

where we introduce F„, Fp and Fe to represent the total momentum transfer rates. Both Fq and F^ are proportional to p and 
so at some value, po, the right hand side of equation 1321 vanishes. Mass elements with densities equal to po will then turn 
around at z and be ejected; if the density is smaller than po then the motion will have been reversed at some height greater 
than z; if it is larger then the mass element will pass by z. Solving equation H32|l for this density gives 



Gn Mbh [FnY + {Fp + Fe){l - Y)] 



(35) 



Clearly po is a function of z and has a maximum value, p*, we name the critical density. Material with a density larger than 
p^, cannot be reversed at any value of z so must be accreted into the bl ack hole. Using Yn = 0.5 these critical densities are 
listed in table (jlj. When a comparison can be made with the results in lFrver fc Meszaroj ll2003l) our results are somewhat 
smaller. 



4 DISCUSSION AND CONCLUSIONS 

The accretion disk model of the central engine of GRBs is an attractive one whose ultimate success will rely on the efficient 
extraction of the energy and momentum of the neutrinos emitted from the disk. The non-spherical geometry of the source 

A least square fit, as in discussed in footnote 1^ gives F^. oc z^-^"^ /{z^ + z'^)^-'^ with = 37.4 km. 
A least square fit gives oc z^-^^ /{z^ + ^'2^-^^ with z* = 40.6 km. 
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makes this easier to achieve than in the context of a spherical source such as in a Type II supernova. The neutrinos are 
essentiahy emitted from a ring: a configuration that allows neutrinos to approach one another at larger interaction angles, 
7, with the concomitant increase in center of mass energy. The bulk of neutrino annihilation occurs in an oblate spheroidal 
region, extending upwards to z ~ 40 km and outward to r ~ 90 km for our model 10, with the highest deposited power in the 
vicinity of the inner edge of the disk. In comparison with other azimuths, the power deposited along the z axis is a minimum. 
The electron/positron pairs would locally thermalize so one should expect a large pressure gradient directed toward the disk 
axis. In addition neutrino absorption and scattering creates a region, close to the black hole, in which the momentum transfer 
is also directed towards the disk axis. Trapping of neutrinos and antineutrinos within the disk leads to significant differences in 
the energy and momentum deposition structure particularly in the region close to the black hole. This may lead to noticeable 
changes in the explosion dynamics. 

These results suggest a scenario whereby the electron/positron pairs created due to neutrino annihilation would, at least 
in part, move initially 'inward and upward' facing little resistance from the e'*'/e~ pairs created along the z axis and accelerated 
by neutrino scattering from behind. At the disk axis the converging e"'"/e~ gas would collide, further increasing the internal 
energy, and then the hot, lepton rich gas flowing rapidly along the z axis. The movement of the e^/e~ gas created by neutrino 
annihilation would have consequences for an explosion. The results in tabled] used only the power deposited along the disk 
axis to calculate the density of material that could be ejected but the movement of the e'^ /e~ pairs from the 'hot spots' 
around the inner disk would suggest that these densities may be too small. The collision at the disk axis of the e^/e~ gas 
from around the torus from would raise the internal energy and, furthermore, the increase in the number density would also 
increase the neutrino opacity leading to a greater transfer of their momentum. Larger critical densities would also imply the 
more prompt formation of a jet. 
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Table 1. Disk Models 



Model 


Disk Calculation 


fa 


a 


a 


Mbh 




an 






[M0S-1] 






[Mq] 


[km] 


[km] 


1 


PWF 


0.01 


0.1 





3 






2 


PWF 


0.01 


0.01 





3 






3 


PWF 


0.01 


0.1 


0.5 


3 






4 


PWF 


0.01 


0.01 





10 






5 


PWF 


0.1 


0.1 





3 






6 


PWF 


0.1 


0.01 





3 






7 


PWF 


0.1 


0.1 


0.5 


3 






8" 


PWF 


0.1 


0.1 


0.95 


6 






gfc 


PWF 


1.0 


0.1 





3 






10 


DPN 


1.0 


0.1 





3 


72.6 


33.6 


11 


DPN 


10.0 


0.1 





3 


240.2 


158.1 



lively. 

It is likely that (unaccounted for) neutrino trapping is important in these models. 



Table 2. Luminosities and Efficiencies 



respec- 



Model 


Ln 

[lO^^ergs-i] 




Lv/r'n 




[10^-"^ ergs- 








1 


0.0094 


5, 


.28 X 10" 


-4 


4.01 X 10" 


-9 


4.25 X 10" 


-7 


2 


0.396 


2, 


.22 X 10" 


-2 


2.90 X 10" 


-6 


7.31 X 10" 


-6 


3 


0.0207 


1, 


.16 X 10" 


-3 


3.24 X 10- 


-8 


1.56 X 10" 


-6 


4 


0.0301 


1, 


.68 X 10" 


■3 


5.28 X 10" 


■9 


1.75 X 10" 


-7 


5 


2.35 


1, 


,31 X 10" 


-2 


4.13 X 10" 


-4 


1.76 X 10" 


-4 


6 


2.14 


1, 


.20 X 10' 


-2 


4.05 X 10- 


-4 


1.90 X 10" 


-4 


7 


5.34 


2, 


.99 X 10" 


-2 


3.58 X 10- 


-3 


6.70 X 10" 


'i 


8 


15.3 


8, 


.57 X 10" 


-2 


4.27 X 10- 


-2 


2.79 X 10" 


-3 


9 


54.3 


3, 


.03 X 10" 


-2 


0.124 




2.29 X 10" 


-3 


10 


109 


6, 


.13 X 10- 


-2 


0.117 




1.07 X 10" 


-3 


11 


1674 


9, 


.37 X 10- 


-2 


4.66 




2.79 X 10- 


3 



Table 3. Peak Momentum transfer per target particle along z axis 



Model ^ Fp ^ Fn/y F^jy Fpiy Fp,y Fgi^ Fe^ 

[gcms~^] [gcms-^] [gcms^'^] [gcms-^] [gcms-^] [gcms-'^] [gcms~^] [gcms^^] 



1 


1, 


.68 


X 


10- 


-17 


1.48 


X 


10- 


-17 


3.71 


X 


10- 


-18 


5.27 


X 


10- 


-18 


3.30 


X 


10" 


-18 


4.69 


X 


10" 


-18 


4.22 


X 


10" 


■19 


2.38 


X 


10- 


■18 


2 


6, 


.15 


X 


10- 


-17 


3.07 


X 


10- 


-17 


1.31 


X 


10- 


-17 


1.08 


X 


10- 


-17 


1.17 


X 


10- 


-17 


9.67 


X 


10- 


-18 


1.54 


X 


10- 


■18 


3.69 


X 


10- 


-18 


3 


9, 


.35 


X 


10" 


-17 


7.54 


X 


10" 


-17 


2.09 


X 


10- 


-17 


2.75 


X 


10- 


-17 


1.86 


X 


10- 


-17 


2.45 


X 


10- 


-17 


2.36 


X 


10- 


■18 


1.37 


X 


10- 


-17 


4 


2, 


.96 


X 


10- 


-19 


2.77 


X 


10" 


-19 


5.93 


X 


10- 


-20 


1.01 


X 


10- 


19 


5.29 


X 


10- 


-20 


9.03 


X 


10- 


-20 


7.43 


X 


10" 


■19 


2.63 


X 


10- 


-20 


5 


1, 


.17 


X 


10- 


-14 


7.11 


X 


10" 


-15 


2.65 


X 


10- 


-15 


2.80 


X 


10- 


-15 


2.35 


X 


10- 


-15 


2.49 


X 


10~ 


-15 


3.02 


X 


10" 


■16 


1.73 


X 


10- 


-15 


6 


3, 


.39 


X 


10- 


-16 


1.27 


X 


10- 


-16 


7.21 


X 


10- 


-17 


4.51 


X 


10- 


-17 


6.41 


X 


10- 


-17 


4.01 


X 


10- 


-17 


8.50 


X 


10- 


■18 


1.57 


X 


10- 


-17 


7 


6, 


.20 


X 


10" 


-14 


3.40 


X 


10" 


-14 


1.40 


X 


10- 


-14 


1.42 


X 


10" 


-14 


1.25 


X 


10" 


-14 


1.26 


X 


10" 


-14 


1.62 


X 


10" 


■15 


9.62 


X 


10- 


-15 


8 


2, 


.97 


X 


10- 


-13 


1.48 


X 


10- 


-13 


6.73 


X 


10- 


-14 


6.79 


X 


10- 


-14 


5.99 


X 


10- 


-14 


6.04 


X 


10- 


-14 


7.96 


X 


10- 


■15 


5.16 


X 


10- 


-14 


9 


1, 


.03 


X 


10" 


-13 


3.03 


X 


10- 


-14 


2.33 


X 


10- 


-14 


1.15 


X 


10" 


-14 


2.07 


X 


10" 


-14 


1.02 


X 


10" 


-14 


2.63 


X 


10- 


■15 


6.41 


X 


10- 


-15 


10 


4, 


.14 


X 


10- 


-13 


1.04 


X 


10- 


-13 


9.30 


X 


10- 


-14 


3.91 


X 


10- 


-14 


8.23 


X 


10- 


-14 


3.47 


X 


10- 


-14 


1.05 


X 


10- 


■14 


2.13 


X 


10- 


-14 


11 


1, 


.72 


X 


10" 


-13 


2.62 


X 


10" 


-12 


3.85 


X 


10" 


-14 


1.03 


X 


10" 


-12 


3.43 


X 


10- 


-14 


9.19 


X 


10- 


-13 


4.38 


X 


10- 


■15 


6.41 


X 


10- 


-13 
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Table 4. Critical Densities; the maximum density of material along the z axis whose motion can be reversed. 



Model 


P* 




L& "-^^^ J 


1 


9.57 X 10"" 


2 


0.215 


3 


1.04 X 10-2 


4 


4.65 X 10-8 


5 


90.8 


6 


4.06 


7 


1.05 X 10^ 


8 


7.26 X 103 


9 


1.51 X 10* 


10 


1.79 X 10* 


11 


8.38 X 105 
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Figure 1. Basis vectors, angles and distances used in tlie calculation of neutrino annihilation, absorption and scattering. The black hole 
is located at the origin and the accretion disk has inner and outer radii of Rin and Rout ■ Neutrinos, emitted from the disk at rn , travel 
along the vector k' toward the point r, z located in the plane formed by i and k. The vector i' lies within the plane containing 6. 

This paper has been typeset from a TJ^X/ M^^X file prepared by the author. 



N 



1 




-1.0x10"^' -5.0x1 0""" 0.0x10"'^ 5.0x10 

r [cm] 



Figure 2. Shows the energy deposited due to neutrino-antineutrino annihilation in units of eVcm~^ s~^. The x-axis is distance radially 
along the disk, while the y axis is vertical distance above the disk. This model is DPN rh = 1 Mq s~^. The density scale height of the disk 
is shown as a solid line while the trapped regions of the neutrinos and antineutrinos are shown as dashed and dotted lines respectively. 
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Figure 3. Shows the momentum transfer from neutrino-electron scattering at all points above the accretion disk in the DPN m = 
1M0 s~^ model. The x and y axis are the same as in Fig. 1^ 
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Figure 4. Shows the momentum transfer from neutrino-proton scattering at all points above the accretion disk in the DPN m = 1 Mq s 
model. The x and y axis are the same as in Fig. |2] 



